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Recent DESY-HERA data on J /<& elasticity distribution show that it emerges mostly as a fast particle. 
Interpreting photoproduction as a collision between a pre-formed charmed hadron and the proton, the out- 
coming J/ty is a leading particle of the collision. We analyse these data using a model formulated to describe 
energy flow in hadron-hadron reactions. The measured J '/vp spectrum can be successfully described in terms 
of this model. We conclude that the observed transparency of the charmed hadron-proton collisions arises 
f> |' because of the particularly small gluonic content of the initial c — c state. 



PACS number(s): 13.85.Qk, 11.55.Jy 



The leading particle effect is one of the most interesting features of multiparticle production in hadron- 
hadron collisions jj]. In these reactions the valence quarks of the projectile emerge from the collision carrying 
the initial state quantum numbers. All produced particles come essentially from the gluons and quark-anti- 
quark pairs already pre-existing in the projectile and target, or radiated during the collision. This qualitative 
picture takes different implementations in the many existing multiparticle production models. In one of them, 
the Interacting Gluon Model (IGM) ^ ||, [|, || , the produced particles (and consequently the energy released in 
the secondaries and lost by the projectiles) come almost entirely from the pre-existing gluons in the incoming 
hadrons. This conjecture may be directly tested using a high energy, nearly gluonless hadronic projectile. In 
this case, according to the IGM, inspite of the high energy involved, the production of secondaries would be 
suppressed (in comparison to the production observed in ordinary hadron induced reactions) and the energy 
would be mostly carried away by the projectile leading particle (LP) which would then be observed with a hard 
xf spectrum. This type of gluonless projectile is available in J/\& photoproduction, where the photon can be 
understood as a virtual c—c pair which reacts with the proton and turns into the finally observed J/\&. There 
are low energy data taken by the FTPS Collaboration |3| and very recently high energy data became available 
at HERA §. 

It has been shown by a number of authors that the simple vector meson dominance mechanism (VDM), 
in which the photon is converted directly to a J /\& before the interaction with the proton, does not describe 
many aspects of data. However, as it was shown in ref. [pi, the charm anti-charm pair can be understood as a 
supperposition of many hadronic states: J/^f, . Whereas relevant for the calculation of the cross sections, 

the inclusion of all the other charmonium states does not change the fact that the projectile, whatever it may 
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be, is a hadron with low gluonic content. For the study of the outgoing J/\& momentum spectrum, this is the 
most important information, which allows us to test the picture proposed by the IGM. 

At lower energies we can compare the momentum spectrum of the J/ty measured in 'yp — > J/ty X collisions 
j(| with the leading meson (jr and K) momentum spectra measured |l| in 7r p — > it X and K p — > K X reactions 
at the same c.m.s. energy. One observes that the charmed leading particles are much harder. This comparison is 
however not completely meaningful because the J/vP data contain a diffractive component which was subtracted 
in the hadronic data. 

At HERA a J/ty z (z — Ej/^/E^) spectrum presumably free from the diffractive component was presented. 
Although the energy is different a comparison with the leading particle spectra measured in hadronic reactions 
is still meaningful because the LP spectra have a relatively weak energy dependence. The J/ VP momentum 
spectrum is clearly harder, according to what we expect in the IGM. The T z spectra, which will be eventually 
measured at HERA in a near future, will be even harder. 

Although there are already some PQCD calculations of J/vP photoproduction j^] the large z values relevant 
here (0.5 < z = Eji^/E^ < 0.9), for which PQCD is not always applicable ||, justify, in our opinion, the more 
phcnomcnological examination of these spectra. In J/\& hadroproduction perturbative calculations fail at large 
Xp, where non-perturbative effects become stronger. Here, these effects may also play some role in the large z 
region. 

We present now our quantitative calculations, which illustrate the qualitative discussion made above. The 
IGM is a model which has been developed primarily to analyse energy-momentum spectra of leading particles 
B 3, ||, ||. This work continues therefore the application of the IGM to photon initiated reactions presented 
in [5] putting them on an equal footing with the hadronic processes studied before |^|, ^ (including diffractive 
dissociation ones 0). 

In Fig. I we show schematically the IGM picture of a photon-proton collision. According to it, during 
the interaction the photon is converted into a hadronic state which interacts with the incoming proton. This 
hadronic state contains the c—c and some gluons and we call it simply "hadron" . The hadron-proton interaction 
follows then the usual IGM picture, namely: the valence quarks fly through essentially undisturbed whereas 
the gluonic clouds of both projectiles interact strongly with each other R In the course of interaction the 
hadronic state looses fraction x of its original momentum and gets excited forming what we call a leading jet 
(LJ) which carries fraction z = 1 — x of the initial momentum. The proton looses fraction y of its momentum 
forming another leading jet. In the IGM we consider two possible types of 7 +p interactions: non-diffractive 
and diffractive. In each of these reactions the J '/\& can come either from the fragmentation of the mesonic 
leading jet or from the hadronization of the central gluonic fireball. In this work we shall concentrate only on 
the data taken at HERA R| for p\ > 1 GeV 2 and 0.5 < z < 0.9. In this case, it is enough to consider the 
single non-diffractively J/*5> produced from the fragmentation of the photonic leading jet, cf. Fig. 1. All other 
contributions can be safely neglected here. 

In order to calculate this spectrum we start (cf., ||, || for details) with the function x{ x i v)> which describes 
the probability to form a central gluonic fireball (CF) carrying momentum fractions x and y of the two colliding 
projectiles: 

x{x,y) : 



2-K^JD xy 
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■ exp \ --L- [{y 2 )(x - (x)f + (x 2 )(y - (y)f - 2(xy)(x - (x))(y - (y))] \ , (1) 



where 



D xy = (x 2 )(y 2 )-(xy) 2 and (x n y m ) = f dxx n f ' dyy m u(x,y). (2) 
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Here xo denotes the normalization factor provided by the requirement that J^dx J^dy x( x > y)0{xy — K^in) = •"■ 
with K m i n = y= being the minimal inelasticity defined by the mass mo of the lightest possible central fireball 
CF (represented by the blob in Fig. 1). The dynamical input of the IGM is contained in the, so called, spectral 
function cu(x,y) given by 

v) = G (y) {*y - K Ln) , (3) 

where G(x) and G(y) denote the effective number of gluons in the charmed hadron and in the proton, respec- 
tively. They are approximated by the respective gluonic structure functions. W is the photon-proton cm. 
energy. a gg is the gluon-gluon cross section, which is computed over a wide range of scales given by xyW 2 . 

x By gluonic clouds we understand a sort of "effective gluonsn" which include also their fluctuations seen as qq sea pairs. 



Whenever the scale is larger than 2.3 GeV lowest order perturbative QCD formulas are used. Otherwise a 
parametrization which represents non-perturbative physics is employed [Q. At the energies W considered here 
the bulk of the interaction happens in the non-perturbative domain, a is the relevant meson-proton cross section. 



The final momentum spectrum of the produced J is then given in terms of x{ x > V) as follows: 



F(z) 



dx / dy x(x,y)S{z 
11 Jo 
1 

dyx(x = 1 - z;y) 

9n.m 



xy 



m^ 



(Mj/v + m ) 2 
W 2 



(4) 



where 



and z 



Vn 



= Max 



ml 



{Mj/y + m ) 5 



(1 - z)W 2 



W 2 
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is the J/^> energy fraction (which in the IGM, where all masses have been consistently neglected 
coincides with the momentum fraction). Because we are dealing here with the leading particle spectra, we 

M x 
and 



have to introduce the additional kinematical constraint, y > ( J - f - , /^+ m °) ; which ensures that the mass Mx 
(Mx = \JyW 1 see Fig. 1) is large enough to produce both the measured J/ty particle of mass Mj/y 
the minimal CF of mass mo, as demanded by the IGM. (In fact, in our case we have to replace Mj/y by 

M T = 

points) 



Mjiy + p^ 7 Px — 1 GeV 2 , to account for the minumum transverse momentum present at our data 



Whenever possible we keep all parameters the same as in the previous applications of the IGM @ - §■ 
However, the inelastic (c—c) — p scattering cross section and the number of gluons in the charmonium state 
are different from the corresponding quantities encountered so far. 

The charmonium-hadron cross section, which we shall approximate by 0j/^_ p , has been subject of intensive 
research in the context of nuclear physics and signatures of quark gluon plasma. Calculations seem to converge 
to <7j™^_ p ~ 6 — 9 mb pfl| . As for the distribution G(x), which we may call G J ^(x), we shall assume that 

it has the same shape as in other mesons, i.e., G J ^(x) — G p ° (x) = G 7r (x) and use, for the latter, the SMRS 
parametrization fTlfl . The specific shape chosen for these distributions does not affect much the results. Their 
normalization p h = J^dxxG h (x), i.e., the amount of momentum allocated to the gluonic component, plays, 
however, a crucial role here and we should try to estimate it somehow. It is known that in a nucleon or in a 
light meson p h ~ 0.5, i.e., gluons carry half of the momentum of the hadron. The charmonium, however, is a 
non-relativistic system and almost all its mass comes from the quark masses. The gluonic field, responsible for 
a weak binding, carries only a small fraction of the energy (and momentum) of the bound state. We expect 
therefore the normalization factor p J /^ of G J ^ (x) to be of the order of the energy stored in the field divided 
by the mass of the state. In the case of a J/ip we have: 



M 



2 m r 



M 



0.033, 



(6) 



where m c has been taken 1.5 GeV and M 



7/V 



3.1 GeV. In the IGM those two parameters are in fact entering 



only as a combination p- ' / '<J % yj% _ p . Because out of those two quantities relevant here only the p J ^ is so far 

' = 9 mb and leave p J ' * to be a free 



completely unknown, we shall, in what follows, take for definiteness C//^_ p 
parameter. 

In Fig. 2 we compare our results (F(z)) with the experimental data. As it can be seen the agreement is very 
good and it was obtained with essentially only one new (but heavily constrained) parameter equal to the ratio 
of the amount of momentum carried by gluons in the cc state and its inelastic cross section with the proton. 
The observed flatteness comes (once o l ^J^_ v is fixed) entirely from the small value of the p' 7 /* parameter, what 

can be seen if we compare our best result (p J ^ — 0.033) shown in the full line, with results for other values 
f p- 7 /*. The dashed line represents the choice p J ^ = 0.066 and the dotted line corresponds to p J /* = 0.016. 
Notice the high sensitivity of the results to the changes of the parameter p" 7 '*. 

To summarize: assuming that the photon may be represented by a hadronic state containing a charm anti- 
charm pair, we have successfully described the leading spectra of photoproduced J/^f's in terms of the IGM. 
At the same time we have demonstrated that (again: within the IGM scheme) they depend crucially on the 
amount of momentum carried by gluons in the charmed hadron, p' 1 ^ (provided its cross section & l J/%_ p with 
the nucleon is known from somewhere else, otherwise they depend of the ratio of these two parameters). It 



means that the knowledge of the leading spectra and the inelastic cross section should allow us to estimate the 
amount of the gluonic momenta in the projectile. For high energies (as those encountered here) this result is 
universal and insensitive to the mass of the projectile under consideration. 
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Figure Captions 

Fig. 1 IGM description of a photon-proton scattering with J/tjj production. 

Fig. 2 Comparison of the IGM distribution F(z) with data of ref. [Q] with restricted acceptance^ > 1 (GeV/c) 2 
and 0.5 < z < 0.9 for fixed value of o~™&_ p = 9 mb and for three different values of 0.066 (dashed 

line), 0.033 (solid line) and 0.016 (dotted line). 
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Figure 2 



